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Abstract 

Despite all the already existing observational data, current models still cannot explain completely the ex- 
cessive energy output and the time variability of GRB. One of the reasons for this is the lack of a good model 
of the central engine of GRB. A major problem in the proposed models with a black hole (BH) in the center is 
that they don't explain the observed evidences of late time activity of the central engine. 

In this paper we are starting the search for a possible model of that central engine as a rotating compact 
body of still unknown nature. The formation of jets in the new model lies entirely on the fundamental Teukolsky 
Master Equation. We demonstrate that this general model can describe the formation of collimated GRB-jets 
of various forms. Some preliminary results are presented. 

1 Introduction 

Gamma-Ray bursts have been mystifying scientists since their discovery in 1969. Their spectral and temporal 
behavior raised many questions offering the unique opportunity to confront our theories with reality. After all 
those years of discoveries and innovations, one thing is clear - GRB cannot be properly understood until there is a 
good model of the central engine. In this paper, we are starting the study of a new model of the central engine of 
GRB that can explain some of their most intriguing characteristics. We do not wont to presuppose a large number 
of hypothetical and unconfirmed by observations specific properties of the central engine. Instead, we prefer to 
focus our consideration on the minimal number of natural assumptions, which seem to be unavoidable in any model 
of central engine, and to investigate their possible consequences. 

1.1 GRB in brief 

Before we present our model, we offer a summary of the properties of GRB that will be important in our studies. 

Gamma-ray burst are explosions on cosmic distances (biggest measured redshift is z > 6.3) that emit huge 
amounts of energy ( ~ 10 51 erg) in very short periods of time (Tgg ~ seconds) ([!]). They have two distinct 
phases - prompt emission and afterglow - produced by different mechanisms. An unexpected feature found by the 
mission SWIFT is the existence of flares ([1], [2] and [3]) - yet not completely explained peaks in the light curve 
superimposed over the continuing afterglow emission. The mechanisms of radiation that are believed to cause the 
burst are synchrotron emission and inverse compton scattering. 

The widely accepted statement us that there exist 2 types of GRB - short GRBs ( Tgo < 2 s) and long GRBs 
(Tgo > 2 s) with some indications of a third type. The temporal separation is supposed to be not only empirical 

- the two types may have different origin and different spectral characteristics. Short GRBs are found in older 
galaxies, they have harder spectrum and they are thought to be a product of binary mergers of compact objects 

- black holes and/or neutron stars. On the other hand, long GRBs are found in young, blue galaxies, with active 
star formation, thus the idea they are a result of the collapse of massive stars([ ]). This hypothesis was confirmed 
by the observation of GRB 060218 which evolved to the supernova SN 2006aj ([4] and [5]). 

Current development of the studies of GRB includes magneto-hydro-dynamical simulations of jets propagating 
in stellar medium with different density or different structure of the jet, testing the type of particles involved, the 
mechanisms of emission and also the type of object that produces them. The observations of GRB continue to 
show interesting characteristics yet to be studied and explained. Because of the dependency of the observations 
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on the characteristics of the detectors, mostly on their energy range, a new space missions are being prepared - 
GLAST 1 among them - that can observe GRB in higher energies which is important for the good resolution of the 
usually hard prompt emission. 



1.2 Fireball model 

Fireball model is one of the most frequently used models in GRB physics [6]. In this model the central engine 
of GRB emits matter in series of shells with different Lorentz factors. When the faster shells catch up with the 
slower ones, the resulting collision called internal shock, produces the hard prompt emission. The deceleration of 
the shells due to the contact with the local medium is believed to produce the softer afterglow. Major set-backs 
of this model are that it cannot explain the late activity of the central engine of which the existence of flares is an 
evidence, nor the nature of those flares. 



2 A Toy Model of Central Engine and Generation of Relativistic Jets 
2.1 General Description 

Our toy model is extremely simple and based on the least possible assumptions. For a start, we have a Kerr black 
hole or some other rotating compact massive object, both described by the Kerr metric ([7]): exactly - in the first 
case, and approximately - in the second case. In Boyer-Lindquist coordinates the metric is: 

ds 2 = (1 - 2Mr/H)dt 2 + AMarsm 2 (6) /T,dtd<f> ~ (E/A)dr 2 - £d<9 2 - sin 2 ((9) [r 2 + a 2 + 2A/a 2 rsin(6»)/S] def 2 (1) 

with A = r 2 - 2Mr+ a 2 , £ = r 2 + a 2 cos 2 6. 

For this metric, Teukolsky studied the linearized perturbations of the Einstein equations and found that the 
equations describing different perturbations generalize to [X]: 
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Following the procedure established by Teukolsky, we perform separation of the variables in the Teukolsky 
equations using the following substitution: 



$ = e^ t+m ^ i S(0)R(r), 



(3) 



where m=...-2, -1, 0, 1, 2... is an integer, lu is the frequency, and S(9) and R(r) are the angular and the radial part 
of the equation. It is important to emphasize that the frequency a; is a complex number: lo = u>r + iu>i. Notice 
that we use different substitution for the time dependence: e lwt , not the original one: e~ luJt , that Teukolsky used. 

One of the most important assumption we use is the stability condition loj > 0, that ensures that the initial 
perturbation won't become infinite with time and it will damp instead. 

Applying the operator (2) L on the function (3) <!> , Teukolsky found that the equations for 8 and r separate 
[8] to an angular and a radial part. The angular Teukolsky Equation (TAE): 
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The Teukolsky radial equation TRE: 

AR lm , rr + 2(s + l)(r - M)R ln>r + V{r)R lm = 0, 
where u = cos(#), A is the constant of separation of the variables and the potential V(r) is: 

uj 2 r 2 (r 2 + a 2 ) +2is (ma (r - M) +uM (r 2 - a 2 )) +2uMr (w a 2 + 2 am) + m 2 a' 2 



(4) 
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x the Gamma-ray Large Area Space Telescope(GLAST ) was renamed to Fermi Gamma-ray Space Telescope after the launch in 2008 
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Here, s=-2, -1, -1/2, 0, 1/2, 1, 2 is the spin weight, related to the value of the spin |s| of the field. It is 
remarkable that all fields can be described with one set of equations. Also, it is important that s Ai m and oj are 
independent parameters. 

To solve that equation, we use the standard notations: r + = M + \jM 2 — a 2 is the event horizon and r_ = 
M — V M 2 — a 2 is the Cauchy horizon. It's easy to see that there is a symmetry between r + and r_ in the TRE 
and that r + and r _ are regular singular points of the TRE, while r = oo is an irregular singular point. 

Using the software package Maple to solve that linear differential equation we acquire two independent exact 
solutions of the radial Teukolsky equation in the outer domain (r > r + ) [9] : 

w(o 2 +r, 2 ) + am ui(a 2 +r 2 ) + am , 

-+i „ / n , r — r. 



R 1 (r)= C t e^ r (r-r + ) r + _r - (r - r_) T +~ r - HeunC [ a,/3^,5,r], „ , „ ) (6) 

and 



R 2 (r) = C 2 e- iur (r-r + y +1 ( r - r_) * ' +1 HeunC (a, -(3, >y,6, r,, + r + ) (7) 

where HeunC is the confluent Heun function (see [10] and [11]) and its parameters are : 
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The angular equation has regular solutions studied by Teukolsky and Press([12]. Beside those solutions, one 
can find polynomial solutions of TAE. As far as we know these polynomial solutions haven't been studied, being 
singular. Explicitly, two solutions of TAE in terms of confluent Heun functions are: 

s£iL(0) = c± ° C ° se (cos(#/2)) |s ~ m| (sin(0/2))- |s+m| x 

/ m 2 - s 2 9\ 
HeunC ±4 fl, \ s—m \, \ s+m \,-4Qs, +2fls~n 2 - A- s,cos - , 



and 

Sj£., m (0) = e ±n c ° s(9) (cos(0/2)) |s - m| (sin(0/2))- |i+m| x 

/ m 2 - s 2 9 
HeunC f =F4^, | s+m |, | s—m\, -4ft s, 2ft s - ft 2 - A - s,sin 2 - 

where f2 — aco. In general, these are two linearly independent solutions of the TAE. They arise from the properties 
of the Heun functions. A special attention is required in the case s £ Z - then s — to 6 Z and the second linearly 
independent solution is not in the above form, but includes an integral of the confluent Heun function. 

2.2 Numerical Results 

Once equipped with the solutions, we proceed with numerical search for the explicit value of the frequency oj. A 
major technical problem in the use of the regular solutions of the angular equation is that we are trying to solve 
connected spectral problem with two complex solutions. It has been solved for the first time by Press & Teukolsky 
([12]) and after that developed trough the mechanism of continued fractions by Leaver [13]. 
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Figure 1: Few examples of jets that we obtained from the solutions of the angular equation. 

In our work, the case is being simplified drastically by the consideration of polynomial solutions, i.e., by the 
assumption that the confluent Heun functions in the angular equation are polynomials. From here, we manage to 
obtain a simple analytic relation A(oj). Using the explicit form of A{ui) we can plot the angular part of the solution 
(3) for certain values of w. 

An important feature of the singular polynomial solutions is that they provide a natural explanation for the 
existence of jets. Examples of such jets in the case s=-l can be seen on figure 1. In that figure, one can see both 
limited and unlimited solutions, although the second one can be constrained adding a form factor with known 
value. Also, we see that different form of the jets arise for different m. Note that the case \s\ — 1, presented in 
this talk, corresponds to perturbations, which describe electromagnetic waves. An analogous treatment of the case 
\s\ = 2, which is more complicated and corresponds to gravitational waves, is presented in [14] and [15]. 

The general solution we use to obtain those plots for \s\ = 1 is: 

S{6) = (cosfl+l) 1 ^ (l-cosfl) 2 ^ e° cos8 HeunC ^40,m+l,m-l,40, ^m 2 - A - 2 fl-fl 2 + i, cos 2 . 

Using the properties of the confluent Heun function and imposing the polynomial conditions on the solutions 
of the angular equation, we obtain the following polynomials, summarized in Table 1. 

The radial equation is much harder to deal with. Using A(ui) and the right boundary conditions for a black 
hole, we want to look for zeroes of the resulting functions. This task is hard considering the profile of the confluent 
Heun functions. 

We use the distributions of the eigenvalues in the complex plane for the singular case: s = — 1, m = 1 for certain 
A(w) (explicitly A s= -i Jm (uj) = —ft 2 — 2 ft m±2 y/QZ + fl m ) to find points resembling to zeroes of the function. After 
that we need an algorithm that can prove they are exactly a zeros of the function and not just a minimum. 

Such new algorithms for searching for the complex roots of a complex function are being developed by the team. 
The results will be reported in the following paper. 

Another problem is that the numerical calculations are additionally burdened by the CPU time required for 
making a complex plot of the solutions by the package Maple. 

3 Conclusion 

Although in very preliminary stage, our model of central engine seems to be able to produce qualitatively some 
of the basic features observed in GRB. The formation of relativistic jets is supposed to be caused by the rotation 
of a compact central body of any nature. At least in a good approximation its exterior gravitational field can be 
described by Kerr solution. The Teukolsky Master Equation (with appropriate boundary conditions) is fundamental 
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Table 1: A table of the functions used to plot Figure 3 



enough to account for both types of GRB with their maybe different origin. The essential assumption we used 
is that the imaginary part of the frequency should be positive. This provides stability of the solutions. Different 
objects can be described by different boundary conditions, though in our case, we used the standard black hole 
boundary conditions for the radial Tcukolsky equation: on the horizon we have only entering waves and on infinity 
we have only outgoing waves. In contrast, we used a novel boundary conditions for the angular Tcukolsky equation 
and this enables us to describe mathematically the collimatcd jets of different forms. Our preliminary results show 
the potential of this mathematical model for description of the central engine as a rotating relativistic compact 
object of any nature. The development of our study will be published in the papers to follow. 
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